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1.  INTRODUCTION 


This  work  is  a  contribution  to  the  solution  of  the  performance  prediction  problem  for 
multiping  detectors  of  target  echoes.  The  formalism  presented  here  applies  to  a  typical  ac¬ 
tive  sonar  detection  problem.  Echoes  resulting  from  pulsed  transmit  waveforms  are  input 
to  a  predetection  stage  (matched  filter) ,  followed  by  a  linear  detector  and  a  postdetection 
integrator. 

Specifically,  we  wish  to  characterize  the  expected  detection  performance  for  an  en¬ 
counter  with  a  real  target,  including  the  dependence  of  target  strength  on  the  observed 
target  aspect  -  that  is,  the  dependence  of  the  signal  amplitude  of  the  received  echo  on  the 
target  aspect.  We  present  an  approach  for  the  generation  of  receiver  operating  character¬ 
istic  (ROC)  curves  for  a  linear  detector  operating  on  signals  in  Gaussian  noise,  where  the 
signal  amplitude  statistics  are  determined  from  empirical  target  strength  data. 

The  results  are  a  practical  extension  of  some  rather  old  results  in  detection  theory. 
Originally,  Marcum^’^  solved  the  detection  problem  for  the  case  of  a  sinusoidal  (constant 
amplitude)  signal  in  Gaussian  noise.  Marcum’s  results  were  extended  by  Swerling,®-^ 
Kaplan,®  Schwartz,®  and  others  to  fiuctuating  signals,  where  the  fluctuations  were  de¬ 
scribed  by  anal3d;ic  densities  that  made  it  possible  to  find  approximate  closed-form  solu¬ 
tions  for  the  multivariate  densities  of  interest.  These  cases  were  reduced  by  Swerling,  in 
work  performed  prior  to  1966  but  only  published  30  years  later,^  to  particular  cases  of 
the  generalized  chi-square  distribution.  In  the  same  work,  he  discussed  the  possibility  of 
using  chi-square  density  functions  to  fit  cases  usually  modeled  with  log-normal  density 
functions. 

These  prior  studies^^^  were  performed  in  the  context  of  radar  physics.  The  empirical 
data  that  motivated  many  of  these  studies  consisted  of  observations  of  target  cross  section 
as  a  function  of  aspect  angle.  In  each  case  the  empirical  probability  density  function  (pdf) 
of  the  measured  data  was  approximated  with  an  analytic  pdf.  This  analytic  pdf  was  then 
used  to  develop  a  solution  to  the  performance  estimation  problem. 

The  methods  developed  in  this  work  use  target  strength  measurements  from  a  target 
of  interest  as  an  empirical  estimate  of  the  signal  pdf  at  the  linear  detector  output.  This 
empirical  pdf  is  then  used  in  the  computation  of  single-ping  and  multiping  detector  ROC 
curves  for  the  target  of  interest.  The  detection  statistic  is  the  envelope  amplitude  (or  sum 
of  envelope  amplitudes)  at  the  linear  detector  output.  The  random  value  taken  by  the 
detection  statistic  depends  on  the  target  aspect  angle,  the  transmit  pulse  duration  and 
bandwidth,  and  the  statistical  characteristics  of  the  additive  noise. 

The  problem  is  to  decide  whether  a  target  is  present  or  absent,  based  on  the  value 
of  a  univariate  random  variable,  the  detection  statistic.  Denoting  the  “target  present” 
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hypothesis  by  and  the  “noise  only”  hypothesis  by  Hq,  the  detection  statistic  z  is 
described  by  the  pdf  pq{z)  if  Hq  is  true,  and  by  pi{z),  if  Hi  is  true.  To  simplify  the 
notation,  the  signal-to-noise  ratio  (SNR)  dependence  is  omitted.  The  detection  statistic 
is  compared  to  a  threshold  value,  b,  that  divides  the  range  of  values  for  i  into  two  regions, 
Rc  =  {z<b}  and  =  {2  >  6}.  The  decision  rule  is  to  decide  that  Ho  is  true  iizeRo,  and 
Hi  is  true  if  2  e  /?i. 

Detector  performance  is  typically  described  by  the  relation  between  the  probability 
of  detection.  Pa,  and  the  probability  of  false  alarm,  for  specified  values  of  the  SNR 
at  the  linear  detector  input.  Pa  represents  the  probability  that  the  detection  statistic  is 
larger  than  a  threshold  value  when  both  signal  and  noise  are  present.  Pja  is  the  probability 
that  the  detection  statistic  is  larger  than  the  threshold  when  only  noise  is  present.  The 
probability  of  detection  is  given  by 


/oo 

Pi{z)  dz  , 


(1.1) 


and  the  probability  of  false  alarm  by 

Pfa  = 

Together  Eqs.  (1.1)  and  (1.2)  represent  the  parametric  form  of  the  receiver  operating 
characteristic  curve.  Implicit  and  explicit  nonparametric  forms,  F{Pfa,Pa)  =  0  and  Pa  = 
Pd{Pfa),  can  be  obtained  by  eliminating  the  threshold  b  between  (1.1)  and  (1.2). 

The  formalism  developed  here  applies  to  both  single-ping  and  multiping  detection 
processes.  In  the  multiping  case,  the  intermediate  observation  vector,  x,  is  made  up 
of  the  linear  detector  outputs  from  successive  echoes,  and  represents  an  observed  value 
of  a  multivariate  random  variable  x.  The  postdetection  integrator  sums  or  averages  the 
detector  outputs  from  successive  echoes  to  generate  a  scalar  quantity,  2,  on  which  the  final 
decision  is  based.  Thus,  the  multivariate  random  variable  x  is  mapped  into  a  univariate 
detection  statistic  2. 

The  empirical  data  used  to  represent  the  signal  amplitude  at  the  linear  detector 
output  consists  of  peak  target  strength  measurements  from  underwater  targets  of  various 
sizes  and  shapes.  The  target  strength  measurements  quantify  the  reflectivity  of  a  particu¬ 
lar  target  as  a  function  of  the  target  aspect  angle  relative  to  the  transmitter  and  receiver. 
These  measurements  were  performed  in  a  controlled  environment  essentially  free  of  noise. 
The  aggregate  noise  at  the  detector  input  (external  reverberation,  receiver  noise,  etc.)  is 
assumed  to  be  zero-mean,  Gaussian,  and  uncorrelated  from  one  echo  to  the  next.  Thus, 
the  noise  amplitude  has  a  Rayleigh  density  at  the  linear  detector  output. 
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The  development  in  the  following  sections  makes  extensive  use  of  the  Rayleigh  and 
Rice  (or  Rayleigh-Rice)  density  functions.  The  Rayleigh  pdf  is  given  by 


and  the  Rice  pdf  is 


p{x,a)  =  ^  exp(-^)  , 


R{r,s,a)  =  ^  exp(-^l^)  /o(^) 


(1.3) 


(1.4) 


The  rest  of  this  report  is  organized  as  follows.  Section  2  develops  the  methodology 
for  computing  ROC  curves  based  on  empirical  target  strength  data  for  the  single-ping 
case.  Although  all  of  the  results  presented  in  Section  2  can  be  obtained  as  particular 
forms  of  the  more  general  multiping  results,  it  is  convenient  to  first  work  out  the  simpler 
single-ping  case.  Section  3  then  derives  the  corresponding  ROC  curve  equations  for  the 
multiping  case.  Section  4  presents  selected  numerical  results  obtained  with  a  MATLAB 
implementation  of  the  analytical  methodology.  These  results  demonstrate  the  capability 
to  provide  an  accurate  quantitative  description  of  performance  for  the  class  of  linear 
detectors  and  for  specific  targets  of  interest. 
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2.  SINGLE-PING  ROC  CURVE  ESTIMATION 
2.1  THE  GENERAL  CASE 

The  target  strength  (signal  amplitude)  is  denoted  in  what  follows  by  s  and  the 
corresponding  random  variable  (RV)  by  s.  The  target  aspect  angle  dependence,  the  source 
of  the  variations  in  the  signal  amplitude,  is  implicit.  The  detection  problem  for  a  given 
aspect  angle  is  equivalent  to  the  detection  of  a  sinusoidal  signal  of  amplitude  s  and  phase  6. 
It  is  assumed  here  that  the  associated  univariates  s  and  6  are  independent  RVs  distributed 
on  [0,oo)  and  [0,27r),  respectively.  The  density  functions  for  s  and  9  are  respectively  /i(r) 
and  fg{6).  The  joint  pdf  for  {§,6)  is  therefore  given  by 

fseM  =  Ms)  f^{e)  .  (2.1) 

If  we  form  the  new  RV 

s  =  (x,y)  =  se*®  =  (s  cos0,  s  sin0)  ,  (2.2) 

the  pdfs  for  (x,  y)  and  (s,  9)  are  related  by 

fxy{x,y)  =  ^fgg{s,9)  =  \fs{s)fe{9)  ,  (2.3) 

s  s 

where  {x,y)  =  (s  cos9,s  sin0). 

Similarly,  the  noise  is  denoted  by  amplitude  n  and  phase  <!>.  The  corresponding  RVs, 
n  and  are  distributed  on  [0,oo)  and  [0,27r)  and  are  assumed  to  be  independent,  with  pdfs 
/ft(n)  and  f^{(t>).  The  joint  pdf  for  (n,^)  is 

Sn^{n,<l>)  =  /ft(n)  f^{4>)  .  (2.4) 

The  new  RV 

{u,v)  —  n  (2.5) 

will  then  have  the  density 

fuv{u,  =  l  4^(n,  4>)  =  ^  fnin)  .  (2.6) 

It  is  assumed  that  s,  0,  n,  and  $  are  mutually  independent  random  variables.  This 
implies  that  s  and  n  are  independent,  and,  consequently,  the  pdf  of  s  +  n  is  given  by  the 
convolution  of  the  pdfs  of  s  and  n. 


fs+n  —  fa*  f A- 


(2.7) 


For  the  single-ping  detection  case,  the  probability  of  detection  and  probability  of 
false  alarm,  (1.1)  and  (1.2),  are  given  by 


/oo 

/ls+n|(r)  dr 


(2.8) 
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and 


Pfaib)  =  fflnlir) 


The  values  of  Pj  and  P/„  in  (2.8)  and  (2.9)  define  a  point  on  the  ROC  curve  for  any 
common  value  of  the  threshold,  b,  corresponding  to  a  specified  SNR  value,  7,  where 


2<t2  • 


(2.10) 


The  problem  is  to  express  /|g+ft|  in  terms  of  the  known  pdfs  fs{s),  f^{e),  fn{n),  and  f^{(j>). 

It  is  assumed  here  that  the  phases  6  and  ^  are  uniformly  distributed  on  [0, 27r)  so  that 


m=fM  = 


(2.11) 


Further,  it  is  assumed  that  n  is  normally  distributed  and  that  each  of  its  components,  u 
and  t),  is  zero-mean  and  has  a  standard  deviation  equal  to  a.  This  implies  that  n  =  |n|  is 
distributed  according  to  the  Rayleigh  pdf 

fn(n)=  p{n,(T)  .  (2.12) 

One  may  note  that  since  n  and  $  are  statistically  independent,  u  and  v  are  uncorrelated 
and  zero-mean.  Given  that  the  RVs  u  and  v  are  Gaussian  and  uncorrelated,  it  follows  that 
they  are  also  statistically  independent. 

The  signal-plus-noise  density  function  of  (2.7)  can  be  written  as  a  convolution  inte¬ 
gral, 

/OO  pOO 

dx'  /  dy'  h{x' ,y')  U{x  -  x' ,y  -  y')  .  (2.13) 

-OO  j  —  OO 

Changing  the  original  integration  variables  to  polar  coordinates  according  to  {x',y')  = 
(r'  cos(9',r'  sinO')  allows  the  right-hand  side  of  (2.13)  to  be  written  in  terms  of  fg  and  /*. 
Similarly,  it  is  convenient  to  express  the  independent  variables  in  polar  coordinates  as 
{x,y)  =  (r  cos0,r  sin0).  With  this  change  of  variables  (2.13)  becomes 


fs+hir  cos9,r  sin0) 


POO  p2Tr 

^  Jo  Jo  (214) 


Using  (2.3)  and  (2.6),  (2.14)  can  be  written  as 


POO  /»27r 

.A(r  cose, r  sin e)  =  dr'  de'fsir')f^{e') 
Jo  Jo 


/n(\/^^  H-  -  2rr'  cos{6'  -  ^)) 
+  r'2  —  2rr'  cos(^'  —  6) 


/.^(tan  ‘( 


1 .  r  sin  g  /  sin  6^  y 
^  r  cos  9  —  r*  cos  0'  ^ 
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(2.16) 


Then  applying  the  uniform  phase  assumption  (2.11)  to  (2.15)  leads  to 


fs+&{r  cos e,  T-  sin  6»)  =  ^  dr'fg{r') 


jn!  /n(Vr^  +  -  2rr-'  cos(0'  -  9)) 

au  - ,  — 


/7>2  _|_  ^/2  _  2rr'  cos(^'  —  6) 


If  the  noise  amplitude  is  described  by  the  Rayleigh  distribution,  the  introduction  of 
(1.3)  into  (2.16)  yields 


1  _ ^ 

fs+nir  cose,r  sin0)  =  2^®  9 


/•°°  .  .  r'2  1  .  T^'co.(e'- 

/  d'r  /i(r')  e-5^  —  /  dO'  e - ^ 

Jo  Jo 


(2.17) 


The  appearance  of  the  Bessel  function 


Io{x)  =  Jo{ix)  =  ^  ^ 


<e~e')d0 


in  (2.17)  makes  it  possible  to  rewrite  this  expression  in  terms  of  the  Rice  pdf,  defined  in 
(1.4),  as 

fs+n{r  cos$,r  sin0)  =  [  ds  f${s)  R{r,s,(7)  .  (2.19) 

27rr  Jq 

Notice  that  fs+&  is  normalized  because  /«  and  R{-,  s,  a)  have  unit  area. 

From  (2.19)  one  can  see  that  the  density  of  s+n  in  polar  coordinates  does  not  depend 
on  the  phase.  This  implies  that  the  RVs  |s  +  n|  and  arg(s  +  n)  are  independent  and  that 

/arg(s+n)(^)  ~  27F’  1^0^ 

fs+A{rcosd,rsme)  =  ^/|s+n|(r)  (2.20) 

and  from  (2.19),  it  follows  that 


noo 

/|8+n|(r)  =  /  ds  fs{s)  Rir,s,a)  . 
Jo 


That  is,  the  unconditional  density  function  f\g+A\  is  the  expected  value  of  a  conditional 
density  function,  namely  the  Rice  function  R{r,s,a),  where  R{r,s,a)  is  dependent  on  the 
constrained  signal  amplitude,  s,  with  univariate  pdf  fs{s). 

One  may  note  that  (2.21)  can  be  expressed  in  terms  of  the  Hankel  transform  defined 

by 

poo 

[Wo/](r)=  /  ds  sf{r)Mrs)  .  (2.22) 

Jo 

Taking  into  account  that  loix)  =  Jo{ix),  (2.21)  can  be  written  as 

/|s+n|(r)  =  ^  ■  (2.23) 

Using  this  formulation,  the  signal-plus-noise  density  may  be  computed  numerically  via 
a  2-D  Fourier  transform,  a  frequency  domain  method  for  computing  the  1-D  Hankel 
transform.*’^® 
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The  probability  of  detection  is  obtained  by  introducing  (2.21)  in  (2.8)  and  evaluating 
the  threshold-dependent  integral, 

roo  poo 

Pd{b)  =  ds  fs{s)  dr  R{r,s,<7)  .  (2.24) 

Since  the  noise  amplitude  is  assumed  to  be  Rayleigh  distributed,  the  probability  of  false 
alarm  can  be  expressed  with  a  simple  analytic  formula.  Introducing  (1.3)  in  (2.9),  the 
probability  of  false  alarm  is  given  by 


Solving  (2.25)  for  b  and  substituting  into  (2.24)  gives 


dr  R{r,s,(T)  . 


(2.25) 

(2.26) 


The  ROC  curve  can  be  represented  graphically  as  a  set  of  points  in  the  plane  with  co¬ 
ordinates  {x,y)  =  {Pfa{b),Pd{b))  by  giving  the  threshold  parameter  b  a  set  of  convenient 
values. 


The  empirical  estimate  of  the  signal  probability  density  function,  fs{s),  is  defined  as 


(2.27) 


where  is  the  set  of  target  strength  observations  at  N  equally  likely  values  of  the 

aspect  angle.  The  empirical  estimate  for  the  probability  of  detection  follows  by  replacing 
the  probability  density  function  in  (2.26)  with  its  empirical  estimate  (2.27).  The  empirical 
estimate  of  Pd  is  then  obtained  as 

1  ^ 

Pd{Pfa,(r)  =  S  PdiPfa,  Sua)  ,  (2.28) 

1=1 

where 

poo 

PdiPfa,s,(T)=  dr  R{r,s,(T)  ,  (2.29) 

In  77^ 

Pfa 

and  where  Pja  is  given  by  (2.25).  Thus,  the  Pd  for  any  given  aspect  angle  can  be  computed 
by  numerically  integrating  (2.29),  and  the  Pd  over  all  aspect  angles  (the  expected  value  of 
Pd  when  the  aspect  angle  is  unknown)  is  computed  from  (2.28). 


2.2  SPECIAL  CASES 

Any  special  case  that  has  an  analytic  solution  is  useful  for  checking  the  numerical 
implementation  of  the  theoretical  results.  Two  such  special  cases,  the  constant  amplitude 
signal  and  the  Rayleigh-fading  signal,  are  considered  below. 
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2.2.1  Constant  Amplitude  Signal 

Expression  (2.21)  takes  a  simpler  form  for  a  constant  amplitude  sinusoidal  signal,  in 
which  case  the  pdf  of  the  signal  amplitude  is  a  Dirac  delta  function, 

fi{r)  =  S{r-ro).  (2.30) 

Prom  (2.21)  it  follows  at  once  that 

/|s+n|(r)  =  R{r,ro,  a)  ,  (2.31) 

which  shows  the  well-known  fact  that  the  pdf  of  the  linear  detector  output  for  a  sinusoidal 
signal  in  Gaussian  noise  is  given  by  the  Rayleigh-Rice  function,  as  defined  by  (1.4).  ROC 
curves  for  this  case  have  been  estimated  and  published  by  Robertson.  “ 

2.2.2  Rayleigh  Distributed  Signal 

Consider  the  case  where  both  the  signal  and  noise  amplitudes  are  Rayleigh  dis¬ 
tributed  at  the  linear  detector  output.  The  signal  and  noise  pdfs  are  given  by 

f${s)  =  p{s,iJ.)  and  U{n)  =  p{n,  a)  ,  (2.32) 

where  p,  and  a  are  the  Rayleigh  parameters. 

Now  (2.21)  can  be  put  in  a  simple  form  with  the  help  of  the  following  relationship 
[9,  eq.  11.4.29,  p.  486]: 

/  =  — 4^0-3^  (3ii/>-l,  5Ja2>0).  (2.33) 

Jo 

Let  1/  =  0,  ,  and  h  =  i^  in  (2.33);  then  one  can  write  (2.21)  as 

=  •  (2.34) 

This  follows  from  the  fact  that  the  sum  of  two  independent  Gaussian  RVs  is  also  a  Gaussian 
RV. 

Substituting  (2.34)  into  (2.8),  it  follows  directly  that  Pd  is  given  by 

ftW  =  exp(4^).  (2.35) 

The  threshold  b  can  be  easily  eliminated  between  (2.25)  and  (2.35)  to  obtain  the  analytic 
relation  between  Pd,  P/a,  and  7  as 

Pd  =  Pfa^  >  (2-36) 

where  7  =  p/a. 
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3.  MULTIPING  ROC  CURVE  ESTIMATION 


3.1  THE  GENERAL  CASE 

We  now  consider  a  processor  that  sums  M  independent  samples  at  the  output  of  the 
linear  detector.  In  particular,  we  consider  the  multiping  processor  where  the  M  samples 
represent  the  peak  linear  detector  output  from  M  successive  echoes.  The  signal  and  noise 
components  at  the  linear  detector  output  for  the  i-th  echo  are 

Hi  =  {nix,niy)  =  (niCos(^i)>^*sin(^i)),  (3.1) 


Si  =  (six,  Siy)  =  (sj  cos(<Ai),  Si  sin((/>i)).  (3.2) 

Each  observation  of  the  multivariate  signal  RV  consists  of  M  scalar  observations  (si, . . . ,  sm)- 
The  ensemble  of  these  observations  form  the  set  {{sik,---,SMk)}k=i-  Introducing  the  nota¬ 
tion 


the  detection  statistic  is  given  by 


Vi  =  Si  +  Hi  , 


M  M 

=  '^Vi  =  ^|Si  -t-Ui 
i=l  i=l 


when  the  target  is  present,  and  by 


M  M 

53ni  =  5]|ni| 

i=l  i=l 


when  only  noise  is  present.  The  Pd  and  Pfa  for  a  given  threshold  b  axe  then 


/OO 

fvi+-+iM{r)  dr 

/OO 

fni-i — yfiM  W  dr  . 


The  SNR  of  an  individual  echo  (ping  i,  aspect  k)  is  defined  by 


^ik 

lik  =  ^  , 


and  the  average  SNR  over  all  aspects  for  a  given  ping  index,  is 


.xzlv' 


Here  we  only  consider  those  detection  scenarios  where  the  noise  statistics  do  not  change 
from  ping  to  ping. 
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As  in  the  1-D  (single-ping)  case,  it  is  assumed  that  the  angular  parts  of  s,  and  n,  are 
independent  of  the  radial  parts,  independent  of  each  other,  and  are  uniformly  distributed 
on  [0,27r).  This  implies  that 

fsii^iiVi)  =  =  —fgf{Si)f^{<j>i)  ,  (3.10) 

oi  Si  ^ 

with  Si  =  |si|  =  y/Six  +  Siy  and  <j>i  =  arctan  and 


(3.11) 


with  m  =  y/Uix  +  Tiiy  and  tpi  =  arctan  The  mutual  independence  of  the  radial  and  angular 
components  and  the  uniform  distribution  of  the  phase  component  on  [0,27r),  for  both  the 
signal  and  noise,  allow  one  to  write 


/siS2,..8Af  (Sl ,  §2, . . . ,  Sa/  )  f (^Ix)  ^ly }  S2xi  ^2yi  •  •  •  >  ^A/xi  ^My) 

M  j 

~  (n  ~)/si’2.-»a;(^1i  ®2)  •••Sm)  ■  02,  •  •  •  ,0A/) 


/  2  \  A/  /  2  \ 

—  y^)  (n  /si.52...«Af  (®1,  ®2,  •  •  •  1  Sa/)  , 


(3.12) 


fniAi-.-Ant  (*^1 ,  *^2,  •  •  • ,  Ua/)  —  /n,  (Ui)  /n2(U2)  •  •  • 


—  /nin2...nA/(ul,n2, .  .  .Ua/)  •  f 

=  — /n,  (ni)/^,,  (0i)  •  •  •  (iPm) 

_  /ni(»i) 

42  27rnj 


Since  the  components  of  n;  are  assumed  to  be  identically-distributed,  independent, 
Gaussian  RVs,  the  pdf  of  nt  =  |ni|  is  Rayleigh  distributed.  The  Rayleigh  parameter  is 
equal  to  the  standard  deviation  of  or  Uiy,  that  is,  cr„,^  =  =  ^i,  and 

fni  =  p{-,(ri),  for  all  i.  (3.14) 

It  follows  therefore  that 

M  ,  y 

/A,n2...AAr(ni,n2,...,nAr)  =  .  (3.15) 

i=l  ^ 

Prom  the  assumptions  it  follows  that  the  multivariates  (si,...,sa/)  and  (ni,...,nA/)  are 
independent.  Consequently,  the  density  function  for  their  sum  is  given  by  the  convolution 
of  their  multivariate  densities. 


/si+fii,...,SAf+fiAf  /si,...,  SAf  *  /ni,...,ftA 
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Using  (3.13),  (3.14),  and  the  notation  Vi  =  Si  +  ni  and  Vi  =  |vi|  =  |si  +  ni|,  an  algebraic 
manipulation  similar  to  that  used  in  the  1-D  case  allows  one  to  express  (3.16)  as 

/v,,...,<rM(vi,...,VM)  =  /  d^S  fs,...Sj^(si,...,SM) 

where  J^Md’^s  =  /o°°  dsi...dsM,  3i+  =  [0,oo),  and  where  3?^  is  the  Mth  Cartesian 

power  of  5t+. 

If  Vi  is  expressed  in  polar  coordinates  as  =  (^;iCosxi,Visinx^),  (3.17)  can  be  written 
as 


where 


/viv.vM  (vi,  •  •  • ,  vm)  =  /vi...vM  (’^I  cosxi,  Wi  sin  Xi,  •  ■  • ,  wm  cosxm,  vm  sinxM) 

~  '  fvi...VM  (^1 1  •  •  •  )  )  fxi--'XM  (xi )  •  •  •  )  Xm)) 

Vi-'-Vm 


f  ^ 

/■Oij.-.jijAf  (^1)  *  ■  *  5  “  /  d  S  7si...sjvf  (^1 » *  ■  *  »  I  I  >  ^3 1  » 

i=l 


/xi-.-XM (xii  •  •  •  > X2)  (27r)^  '  (3.20) 

Thus,  under  the  specified  assumptions,  the  radial  and  angular  components  of  (vi,...,vm) 
are  mutually  independent  and  each  angular  component  is  uniformly  distributed  on  [0,27r). 

The  density  of  a  sum  of  univariates  z  =  vi+v2-\ - \-vm  is  determined  by  their  joint 

density  according  to 

/^(^)  =  /  d^~^W  fi^e2...VM{‘Wl>'^2,---,'WM-l,Z-Wi-W2 - WM-l)  ■  (3.21) 

The  expression  (3.21)  is  the  generalization  of  the  well-known  result^®  that,  given  two  RVs 
X  and  y  with  joint  pdf  f&y{x,  y),  the  pdf  of  their  sum,  z  =  x  +  y,  is  given  by 


/OO 

U^{z-w,w)dw.  (3.22) 

-OO 

Thus  (3.21)  becomes  the  M-fold  convolution  of  the  1-D  marginal  densities  of  the  univariate 
RVs  when  the  individual  RVs  are  independent. 

If  (3.19)  is  introduced  in  (3.21)  and  the  order  of  integration  is  changed,  the  M- 
dimensional  integration  in  (3.21)  will  involve  only  the  product  of  the  Rice  densities  from 
(3.19),  and  will  represent  the  convolution  of  these  Rice  functions.  Formally,  the  result  is 

Mz)=f  d^s  fs^...s^{si,...,SM)  G(z,si,...,SM,cr)  ,  (3.23) 
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where 


G(2,Si,...,s„,ct)  =  ,  (3.24) 

and  where  represents  the  M-fold  convolution  operator.  Similarly  to  (2.21)  for  the  1-D 
case,  (3.23)  shows  that  the  unconditional  pdf,  fi{z),  is  the  expected  value  of  a  conditional 
density  function,  namely  the  convolution  of  Rice  functions  given  by  (3.24).  The  conditional 
density  function  (3.24)  is  dependent  on  the  constrained  signal  amplitudes,  (si, . . .,sm),  that 
are  distributed  according  to  the  multivariate  pdf  /i(s). 

The  probability  of  detection  follows  from  (3.6)  and  (3.23)  as 

Pd{b)  =  P^{b,Su...,SM,(T)  ,  (3.25) 

where 

r<X) 

Pd\b,si,...,8M,<^)=  I  dz  G{z,Si,...,SM,a)  .  (3.26) 

*fb 

The  expression  (3.25)  shows  that,  for  a  specified  threshold,  the  probability  of  detecting  a 
fluctuating  target  by  summing  M  echos  is  the  expected  value  of  the  constrained  probability 
of  detection  of  M  “point  targets”  with  the  corresponding  ensemble  of  “constant”  target 
strengths  distributed  according  to  the  multivariate  pdf 

The  probability  of  false  alarm  can  be  derived  from  (3.25)  by  setting  the  signal  com¬ 
ponents  to  zero.  Thus  the  signal  pdf  is  given  by 

M 

fsi...SM{si,...,SM)  =  YlS{si)  ,  (3.27) 

i=l 

and  this  leads  to 


fOO 

Pfa{b)  =  Prf  (6,  Si  =  0, . . . ,  SAr  =  0,  <r)  =  /  dz  cr)](2:)  ,  (3.28) 

J  f) 

where  p{r,a)  is  the  Rayleigh  density  (1.3).  This  could  have  been  written  directly,  since 
the  noise  is  assumed  to  be  Rayleigh  distributed  and  independent  between  echos.  Together 
(3.25)  and  (3.28)  constitute  a  parametric  expression  of  the  ROC  curve. 

The  empirical  estimate  for  the  probability  of  detection  follows  from  (3.25)  if  the  pdf 
/*,...»M(si,-..,SAf)  is  replaced  by  the  empirical  multivariate  pdf 

^  N  M 

/iT..»A,(si.--MSM)  =  ^^JJ(5(si-Sifc)  .  (3.29) 

^  k=l i=l 

The  result  is 

Pd{b)  =  j;rYl  /  =  —  J]p«(6,siA,...,SA/fc,<T)  .  (3.30) 

fc=i  •'  fc=i 
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3.2  SPECIAL  CASES 


Formally  simpler  results  can  be  obtained  if  the  pdf  for  the  sequence  of  signal  am¬ 
plitudes  has  a  manageable  analytical  form.  Such  is  the  situation,  for  example,  when 
the  signal  amplitudes  are  independent  from  ping-to-ping  (random  target  aspect),  when 
the  signal  amplitudes  are  identical  ping-to-ping  (no  aspect  change),  or  when  the  signal 
amplitude  is  constant. 

3.2.1  Independent  Signal  Amplitudes  Ping-to-Ping 

Given  empirical  target  strength  data  as  a  function  of  aspect,  the  independent  echo 
case  would  correspond  to  a  random  change  in  aspect  between  echos.  While  this  case  is 
nonphysical,  it  is  useful  in  checking  the  code  developed  to  estimate  the  multiping  ROC 
cmwes. 

If  the  echos  are  independent,  then  the  pdf  for  the  ensemble  of  amplitudes  is  the 
product  of  the  marginal  pdfs, 

/I1...SM  (^1)  •  •  •  )  ®A^)  ~  (®l)  /s2(®2)  ■  ■  ■  /sM(^Ar)  )  (3.31) 

and  (3.19)  can  be  written  as 

^  1*00 

fvi . =  n  /  dSi  fs^isi)  R{vi,Si,cr)  .  (3.32) 

i=i  Jo 

If  this  expression  for  is  introduced  in  (3.21),  it  follows  that  the  pdf  for 

the  sum  of  amplitudes  is  the  convolution  of  M  identical  replicas  of  the  univariate  pdf  for 
the  amplitude  of  one  echo  given  by  (2.21),  i.e., 

Mz)  =  [*,^l/|s+n|](^)  .  (3-33) 

This  expression  could  have  been  written  immediately,  based  on  the  independence  of  the 
signal  and  noise  components. 

The  empirical  estimate  of  fziz)  follows  from  (3.46)  if  /|s+i>|(2)  is  replaced  by  its  em¬ 
pirical  estimate 

/|f+n|(^)  = 

i=l 

obtained  from  (2.21)  and  (2.27).  The  empirical  probability  of  detection  can  now  be  esti¬ 
mated  numerically  by  using  (3.34),  (3.33),  and  (3.6). 

3.2.2  Identical  Signal  Amplitudes  Ping-to-Ping 

If  the  observed  target  strength  (signal  amplitude)  is  constant  for  any  set  of  M  suc¬ 
cessive  pings,  then  the  pdf  of  the  M- variate  signal  is 

/si...sm(s1)-  •  •  ,Sm)  =  /li(si)  S{S2  -  Sj)  •  •  -5(5^1/  “  Si)  •  (3.35) 
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If  (3.35)  is  used  in  (3.19),  the  result  is 


fviV2  --^M  (^1 


f>oo  M 

,V2,.  ..,Vm)  —  ds/s(s)p  J?(i;j,s,cr)  , 


where  R  is  the  Rice  function  (1.4)  and  si  is  simply  denoted  by  s.  If  (3.36)  is  introduced 
in  (3.21)  and  the  order  of  integration  is  changed,  the  result  is 


fzi^)  =  f  ds  fs{s)  G(z,s,<t)  , 
Jo 


where  G{z,s,a)  is  the  convolution  of  P  identical  replicas  of  the  Rayleigh-Rice  pdf  with 
respect  to  the  first  dependent  variable,  i.e., 

G{z,s,cr)  =  (*l^ii?(-,s,a)](z)  .  (3.38) 


Recalling  that  z  —  vi  +  +  •••  +  v^,  the  probability  of  detection  follows  from  (3.6) 

with  (3.37)  and  (3.38)  as 


/OO  poo 

dr  J  ds  fs{s)  G(r,s,a) 

poo  pOO 

=  /  ds  fs{s)  /  dr  G{r,s,a) 
Jo  Jb 

poo 

=  ds  h{s)P^{b,s,a)  . 

Jo 


The  expression  P^(b,s,<T)  =  dr  G{r,s,a)  is  the  single-ping  version  of  (3.26)  and  has  been 
tabulated  by  Robertson.”  This  special  case  corresponds  to  a  random  (unknown)  aspect 
target,  but  where  the  aspect  remains  fixed  over  any  given  set  of  M  pings  (observations). 
Expression  (3.39)  shows  that  the  probability  of  detection  (.Pd)  for  this  special  case  is  given 
by  a  weighted  average  of  the  Pa  for  M  constant  amplitude  signals  where  the  weighting  is 
defined  by  the  pdf  of  the  signal  amplitudes,  /«(«). 

The  empirical  estimate  for  the  probability  of  detection  is  obtained  by  replacing  the 
pdf  of  s  in  (3.39)  with  its  empirical  estimate. 


ffis)  =  ^f2^is-s,). 


(3.40) 


Here  {sfc}£Li  corresponds  to  the  set  of  target  strength  values  (signal  amplitudes)  measured 
for  a  uniformly  distributed  (equally  likely)  set  of  aspect  angles.  The  result  is 


Pdb)  =  ^'ZPhb,Sk,a)  . 


(3.41) 
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3.2.3  Constant  Amplitude  Signals 

If  the  signal  amplitude  for  each  ping,  At,  is  a  constant,  A,  then 

=  Sisi  -  A)---5{sm  -  A)  .  (3.42) 

With  (3.42),  (3.19)  becomes 

M 

(ui,...,vm)  =  JJii(t;i,A,(T)  ,  (3.43) 

i=l 

and  (3.21)  reduces  to 

f,{z)  =  [*fi,R{;A,<T)]iz)  .  (3.44) 

This  could  have  been  obtained  from  (3.38)  by  making 

fsis)  =  6{s-A).  (3.45) 

Substituting  (3.44)  into  (3.6),  the  resulting  probability  of  detection  is 

poo 

Pl{b)  =  l  dz  [*fiiRi;A,<7)]{z)  .  (3.46) 

The  ROC  curves  derived  from  (3.46)  and  the  corresponding  Pfa  (3.28)  were  charted  by 
Robertson. 
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4.  NUMERICAL  EXAMPLES 


This  section  provides  selected  examples  of  ROC  computations  to  illustrate  the  ana¬ 
lytical  development  presented  in  the  preceding  sections.  Equations  (3.25)  and  (3.28)  allow 
straightforward  numerical  computation  of  linear  detector  ROC  curves  for  analytical  pdfs 
as  well  as  for  empirical  data. 

For  the  case  of  a  constant  target  (constant  signal  amplitude)  and  independent,  Gaus¬ 
sian  noise,  the  Pj.  can  be  computed  by  convolving  the  Rician  pdf  for  the  desired  SNR  with 
itself  M-1  times,  where  M  is  the  number  of  pings  to  be  integrated  [see  (3.46)].  The  cor¬ 
responding  Pfa  is  given  by  (3.28).  Example  ROC  curves  for  the  constant  signal  case  are 
shown  in  Figure  4.1  for  SNR  =  9  dB  and  for  the  integration  of  1,  2,  and  4  echoes. 

The  numerical  computations  for  the  case  of  a  constant  amplitude  target  were  vali¬ 
dated  through  comparisons  with  equivalent  results  published  by  Robertson^^'^^  for  selected 
test  cases.  The  Pd  values,  for  a  given  Pfa  and  SNR,  computed  using  the  methods  devel¬ 
oped  here,  differed  from  the  corresponding  Pd  values  given  by  Robertson  by  less  than  1%. 
Equivalently,  the  small  differences  observed  in  the  Pd  values  correspond  to  a  difference  of 
less  than  0.05  dB  in  the  SNR  required  to  achieve  a  specified  Pd  and  Pfa- 

ROC  curves  for  a  Rayleigh-fluctuating  target  in  Gaussian  noise  are  shown  in  Figure 
4.2.  The  signal  amplitudes  are  assumed  to  be  independent  from  ping-to-ping.  Curves 
are  shown  for  1-,  2-,  and  4-ping  integration  and  an  average  SNR  of  9  dB.  The  average 
SNR  is  defined  as  a^f2a^,  where  a  is  the  Rayleigh  distribution  parameter  and  a  is  the 
standard  deviation  of  the  Gaussian  noise.  The  pdf  of  the  multiping  detection  statistic  in 
(3.23)  is  obtained  by  computing  the  normalized  sum  of  Rayleigh-weighted  Ricians  for  the 
single-ping  case  (2.21)  and  then  convolving  the  resulting  pdf  with  itself  the  appropriate 
number  of  times. 

The  numerical  computations  for  the  Rayleigh-fluctuating  target  ROC  curves  were 
validated  through  comparison  with  the  closed-form  solution  for  the  single-ping  case,  eqs. 
(2.35)  and  (2.36),  and  with  the  results  of  a  Monte  Carlo  simulation  for  the  multiping  case. 
The  reduction  in  performance  (lower  Pd  for  a  given  Pfa)  for  a  Rayleigh-fluctuating  target 
with  average  SNR  =  9  dB,  relative  to  the  performance  for  a  constant  target  with  SNR  = 
9  dB,  is  evident  from  comparison  of  Figure  4.2  with  Figure  4.1. 
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Figure  4.1  Example  ROC  curves  for  a  constant  amplitude  target. 
SNR  =  9  dB;  1-,  2-,  4-ping  integration 


Figure  4.2  Example  ROC  curves  for  a  Rayleigh-fluctuating  target. 
Avg.  SNR  =  9  dB;  1-,  2-,  4-ping  integration 
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Now  consider  the  estimation  of  detector  performance  for  two  actual  targets,  referred 
to  here  as  target  A  and  target  B.  Figures  4.3  and  4.4  compare  the  measured  target  strength 
values  for  the  two  targets.  Figure  4.3  shows  the  polar  scattering  patterns  (relative  target 
strength  as  a  function  of  target  aspect  angle),  while  Fig.  4.4  shows  the  cumulative  distri¬ 
bution  functions  (cdfs)  of  the  target  strength  values  over  all  aspect  angles.  As  seen  from 
Fig.  4.3,  the  target  strength  of  target  A  is  larger  than  that  of  target  B  for  most  aspect 
angles.  In  Fig.  4.4,  the  cdf  for  target  B  has  been  shifted  by  approximately  -1-7  dB  relative 
to  the  cdf  for  target  A  in  order  to  make  the  average  SNR  values  equal  for  the  two  targets. 
The  target  strength  characteristics  for  these  two  targets  are  clearly  very  different. 

Figure  4.5  illustrates  how  the  expected  detection  performance  over  all  aspects,  i.e., 
the  expected  detection  performance  for  a  random  aspect  encounter  with  the  target,  de¬ 
pends  on  the  target  strength  characteristics  of  these  two  targets.  ROC  curves,  computed 
using  (3.28)  and  (3.30),  are  shown  for  both  single-ping  detection  and  for  the  integration  of 
four  pings.  The  target  strength  values  and  noise  statistics  were  scaled  so  that  the  average 
SNR  was  9  dB  for  both  targets,  where  the  average  SNR  is  defined  by  Eq.  (3.9).  Thus, 
the  differences  between  the  ROC  curves  for  targets  A  and  B  result  from  the  differences  in 
the  shapes  of  the  target  strength  cdfs  (Fig.  4.4)  rather  than  from  absolute  differences  in 
the  target  strength  values  (Fig.  4.3)  for  these  two  targets. 

In  many  applications,  detection  performance  against  a  particular  target  is  computed 
based  on  a  point  target  representation;  i.e.,  a  single  target  strength  value  is  used  to 
represent  the  (statistical)  encounter  with  the  target.  A  common  approach  is  to  use  an  order 
statistic,  from  the  target  strength  versus  aspect  values,  that  corresponds  to  the  desired 
probability  of  detection.  For  example,  if  a  90%  probability  of  detection  is  desired,  the 
lower  10th  percentile  target  strength  value  is  used  to  represent  the  target.  The  objective 
of  this  procedure  is  to  force  a  conservative  estimate  of  performance  prediction  relative  to 
the  actual  performance. 

Figure  4.6  compares  the  ROC  curve  for  target  B,  computed  using  the  methods  of 
this  paper  that  incorporate  the  target  strength  measurements  for  all  aspects,  with  ROC 
curves  for  a  point  target  (constant  signal  amplitude)  assumption.  The  blue  curve  shows 
the  estimated  ROC  curve  for  target  B  for  4-ping  integration  and  an  average  SNR  of  10.8 
dB.  The  red  curves  are  the  corresponding  4-ping  ROC  curves  for  the  constant  (point) 
target  assumption  at  SNRs  equivalent  to  the  10th  (lower  curve)  and  50th  (upper  curve) 
percentile  values  of  the  target  strength  distribution  for  target  B.  The  lower  10th  percentile 
corresponds  to  an  SNR  of  2.5  dB  while  the  50th  percentile  SNR  is  7.0  dB.  Clearly,  the 
ROC  curves  for  the  point  target  assumption  provide  a  poor  representation  of  the  expected 
detection  performance  for  this  target. 
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Cumulative  Distribution 


Figure  4.3  Relative  target  strength  (dB)  versus  aspect 
for  two  mine  targets;  target  A  (red),  target  B  (blue) 


Relative  Target  Strength  (dB) 

Figure  4.4  Cumulative  distribution  of  target  strength  values; 
target  A  (red),  target  B  (blue) 
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Probability  of  False  Alarm 

Figure  4.5  ROC  curve  comparison  for  targets  A  (red)  and  B  (blue) 
Avg.  SNR  =  9  dB;  4-ping  integration 


Probability  of  False  Alarm 


Figure  4.6  ROC  curve  comparison  for  target  B  (blue)  with 
constant-signal  point  targets  (red) 


Figure  4.7  shows  a  set  of  ROC  curves  for  target  B  for  the  case  of  2-ping  integration 
for  average  SNRs  of  6,  9,  12,  15,  and  18  dB.  Figure  4.8  shows  a  similar  set  of  ROC  curves 
for  target  B  with  an  average  SNR  of  9  dB  where  the  number  of  pings  integrated  varies 
from  1  to  16  pings.  While  these  curves  assume  a  random  target  aspect  for  the  encounter 
(all  aspects  are  equally  likely),  it  is  straightforward  to  include  any  a  priori  knowledge 
about  the  target  aspect  or  range  of  target  aspects  in  the  ROC  curve  computations. 

The  ROC  curves  in  the  previous  figures  show  the  expected  multiping  detection  per¬ 
formance  when  all  target  aspects  are  equally  likely  and  when  the  aspect  is  constant  over 
any  given  set  of  multiping  observations.  The  effect  of  ping-to-ping  changes  in  aspect  on 
detection  performance  is  illustrated  in  Fig.  4.9.  ROC  curves  are  shown  for  target  B  for 
an  average  SNR  of  9  dB  with  8  pings  integrated  for  a  fixed  aspect  (no  aspect  change),  a 
shift  in  the  observed  aspect  of  0.4  degrees  per  ping,  and  a  shift  of  2.0  degrees  per  ping. 
In  addition,  the  ROC  curve  is  shown  for  the  case  of  a  random  aspect  shift  from  ping  to 
ping.  As  in  the  previous  figures,  the  ROC  curves  represent  the  expected  performance  for 
a  random  aspect  encounter  with  the  target.  The  independent  aspect  case  represents  an 
approximate  upper  bound  on  performance  for  ping-to-ping  aspect  changes  typical  of  a 
real-world  encounter. 


Probability  of  False  Alarm 


Figure  4.7  Target  B  ROC  curves  versus  SNR. 
2-ping  integration;  Avg.  SNR  =  6,  9,  12,  15,  18  dB 
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5.  SUMMARY 


This  work  presents  a  method  for  incorporating  target  scattering  patterns  directly 
into  the  computation  of  statistical  detection  performance  for  a  linear  detector  followed 
by  a  post-detector  integrator  that  sums  the  detected  peaks  over  multiple  echoes.  Tar¬ 
get  scattering  patterns  (e.g.,  acoustic  target  strength  or  radar  cross-section)  characterize 
the  echo-to-echo  fluctuations  in  level  due  to  changes  in  observed  target  aspect.  In  this 
work,  the  scattering  patterns  of  interest  are  obtained  from  direct  measurement  of  acoustic 
target  strength  of  underwater  bottom  mines;  however,  scattering  patterns  obtained  from 
theoretical  models  or  other  means  can  be  used  as  well. 

As  illustrated  by  the  examples  in  the  previous  section,  the  methodology  presented 
here  supports  estimation  of  multi-ping  Receiver  Operating  Characteristic  (ROC)  curve 
performance  for  any  simulation  scenario  of  interest.  For  example,  detection  performance 
can  be  estimated  for  a  random  encounter  with  a  specific  target  (i.e.,  initial  aspect  un¬ 
known),  for  a  specific  sonar-target  encounter  geometry  as  a  function  of  the  initial  observed 
target  aspect,  or  as  a  function  of  the  range  of  target  aspects  observed  in  a  sonar-target 
encounter.  Alternatively,  the  method  can  be  used  when  the  received  signal  amplitudes  are 
described  by  an  analytic  distribution.  The  methodology  can  be  extended  to  the  general 
detector-averager  problem  for  sampled  time-varying  echo  envelopes  given  that  successive 
samples  of  the  received  signal-plus-noise  are  independent. 
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